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I. Introduction 
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One of the most difficult problems the manufacturing community has faced during recent years 
has been to accurately assess the physical state of anisotropic high-performance materials by 
nondestructive means. Destructive testing of materials, although it supplies useful information for 
the design engineers in the development stage of a material, is inappropriate for ascertaining the 
state of the material during assembly testing and routine maintenance procedures. Conventional 
nondestructive ultrasonic measurement techniques, currently employed with much success, can be 
improved to yield even more information about the physical state of a material. In the past, 
communication between the physical science and engineering communities has been somewhat 
inadequate. Each community has valuable information to contribute towards the assessment of 
material integrity. Measured ultrasonic parameters can be related to the common engineering 
parameters to yield information useful to both communities. 


‘ In order to advance the design of ultrasonic nondestructive testing systems, a more 
fundamental understanding of how ultrasonic waves travel and interact within the anisotropic 
material is needed. The relationship between the ultrasonic and engineering parameters needs to be 
explored to understand their mutual dependence. One common denominator is provided by the 
elastic constants. Accurate measurements of these physical parameters can yield information useful 
in determining the physical integrity of the material. With this in mind, advanced ultrasonic 
measurement systems can be designed with the internal nature of the material in mind. The type of 
material and the very question the investigator is trying to answer should play a major role when 
setting up a measurement system. 

In Section II of this Progress Report we discuss our preparation of specific graphite/epoxy 
samples to be used in the experimental interrogation of the anisotropic properties (through the 
measurement of the elastic stiffness constants). Accurate measurements of these constants will 
depend upon knowledge of refraction effects as well as the direction of group velocity propagation. 
Section III discusses our continuing effort for the development of improved visualization 
techniques for physical parameters. Group velocity images are presented and discussed. In order 
to fully understand the relationship between the ultrasonic and the common engineering 
parameters. Section IV discusses the physical interpretation of the linear elastic coefficients (the 
quantities that relate applied stresses to resulting strains). This discussion builds a more intuitional 
understanding of how the ultrasonic parameters are related to the traditional engineering 

parameters. 
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II. Sample Preparation and Preliminary Velocity Measurements 

As discussed in previous Progress Reports, uniaxial graphite/epoxy laminates exhibit 
hexagonal symmetry. Five distinct elastic stiffness coefficients are required to describe the elastic 
properties of this material. Four of the five constants may be measured by propagating 
longitudinal and shear waves perpendicular and parallel to the fiber orientation. The fifth elastic 
constant must be measured at an angle that is neither parallel nor perpendicular to the fiber 
direction. In our Progress Report of 9/86 a discussion of the existence of a possible optimum 
angle at 77° with respect to the fiber orientation, for measurement of this fifth elastic constant, was 
predicted based on analysis of the propagation of errors. 

We have prepared three uniaxial samples, received from NASA Langley Research Center, for 
the measurements. Two of the samples were surface ground so that their sides were parallel and 
perpendicular to the fiber orientation as illustrated in Figure 1. 



Figure 1 : Two samples prepared for the propagation of longitudinal and shear 
waves parallel and perpendicular to the fiber orientation. 

The third sample was prepared so that insonification normal to the surface will produce ultrasonic 
waves whose phase velocity will propagate inside the sample at an angle of IT with respect to the 
fiber orientation. Figure 2 illustrates the two steel stages machined and used to prepare the third 
sample of graphite/epoxy. Also illustrated in Figure 2 is the finished sample. 

Time resolution of the received ultrasonic signals requires an understanding of refraction 
effects and possible mode conversion at interfaces. An understanding of the directions of phase 
and group velocity propagation inside the material is also required (see Section III of this Progress 
Report). Using results obtained from linear elastic theory we are evaluating a preliminary design 
of the measurement system. A full report of these investigations will be forth coming. 
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Stage 1 


Stage 2 


Side View 




Final Sample Side View 



Figure 2: Preparation of a uniaxial graphite/epoxy sample for measurement of 
the fifth elastic constant 
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III. Visualization of Group Velocity Surfaces for Uniaxial Graphite/Epoxy 
Composites 

For anisotropic materials the phase and group velocity need not be collinear. The group 
velocity is defined in terms of the modulation of the wavefronts. Consider an ultrasonic wave 
emitted from a conventional planar cylindrical cross section transducer. The wavepacket is 
restricted in two dimensions by the diameter of the transmitting transducer and in the third 
dimension by the ultrasonic pulse length. The wavefronts propagate in the direction of k (along 
the axis of the transducer) while the wavepacket propagates along the direction of group velocity 
(modulation envelope). If the phase and group velocities are not collinear, the receiving transducer 
will have to be offset with respect to the phase velocity direction in order to intercept the ultrasonic 
pulse. By using the information obtained from linear elastic theory we are able to predict the group 
velocity direction for a given initial phase velocity direction. This information will enable the 
investigator to incorporate such knowledge into the design of the measurement system. Although 
numerical calculations produce useful numbers, they do very little to raise the physical intuition of 
the investigator about the intrinsic nature of the material with which he is working. Our hypothesis 
is that a three dimensional representation of the group velocity surfaces for each allowed ultrasonic 
mode of propagation would greatly improve the insight one has when designing an advanced 

measurement system. 

Using the fact that the fiber direction is an axis of six-fold symmetry, which implies that 
orthogonal planes to this axis are transversely isotropic, we have produced three dimensional 
group velocity surfaces for uniaxial graphite/epoxy composites. Figure 3 displays the group 
velocity surfaces for the three propagation modes. 

The way to interpret these surfaces, for each propagation mode, is to place the observer inside 
the surface at the origin. When he/she looks in any direction, the length of the vector from the 
origin to a point on the surface represents the magnitude for the group velocity that propagates in 
that direction. In general for anisotropic materials, this does not necessarily result from the 
propagation of the phase velocity in this direction. We are currently working on a representation 
which will enable the viewer to more easily visualize the relationship between the initial phase 
velocity direction (i.e., the direction in which the wave is launched) and the resulting group 
velocity direction. 



- 5 - 



- 6 - 


IV. Physical Interpretation of the Elastic Matrices 

In the previous Sections of this Progress Report we have stressed the importance of 
understanding the intrinsic nature of the material that is to be evaluated. The more we understand 
about the intrinsic physical properties of a material, the more intelligent the design of the 
measurement system can be. The ultimate goal is the assessment of material integrity by a 
nondestructive means. Ultrasonic measurement techniques have proven to be a very attractive 
approach for solving this problem. Information about the strengths and weaknesses of a material 
can be derived from ultrasonic parameters with the added feature that these parameters can be 
directly related to the traditional engineering parameters, such as Young's moduli, Poisson's ratios 

and the like. 

In this Section the relationship between the stress and strain in a linear elastic material, via the 
elastic stiffness and compliance matrices, is discussed. A physical interpretation of the form of the 
elastic matrix as well as the individual components is presented. The relationship between the 
elastic coefficients and the Young's moduli, Poisson's ratios, and the pure-shear moduli is also 

discussed. 

The Constitutive Relation: Hooke's Law 

The constitutive relation, relating stress and strain, in the most general case requires 81 
coefficients to describe the elastic properties of a linear elastic solid. For a lossless material the 
stress (Ty) is proportional to the strain (Ski) via the real-valued elastic stiffness coefficients (Cijki)- 

3 3 

T ij = X X Ci i kl Skl (i,j = 1,23) 

k = t 1=1 

or T = c : S (1) 

Stress denotes a force acting on an area. It is useful to consider the stress tensor in dyadic form to 
explicitly show the direction of the applied force and the surface area to which it is applied. This 
representation, common in engineering literature, yields a more physical interpretation of the stress 
tensor. A dyad is defined as the direct multiplication of two vector quantities to produce a nine 
component 2nd-rank tensor quantity. We can write the stress tensor in dyadic form as 

T = = <F,.F 2 .F 3 ><i -fc.J? (a direct product) 


or in matrix form as 
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[Tijl 


Fi_ Fj_ Fi_ 

A} A 2 A 3 

f 2 F 2 F2_ 

“ A 3 where F; is the component of force in the ith direction 

F 3 F 3 F 3 and Aj is the area of the surface whose normal is in the 

Ai A 2 A 3 J jth direction . (3) 


In Figure (4) the notation for the applied stresses on a unit cube is illustrated. The first 
subscript denotes the direction the force is acting; the second denotes the direction of the normal 
vector to the face of the cube on which the force is applied. 



Figure 4: The stress components acting on each face of a unit cube within a stressed 
body (ijkl system). 

Figure 4 along with the previous discussion demonstrates the necessity for the stress being a 
double subscript quantity. 

From classical linear elastic theory, the maximum number of distinct elastic constants for a 
material is reduced from 81 to 36 due to the symmetry of the stress and strain tensors. That is, 

Tij = Tji and S y = % 
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from which it necessarily follows that 

Cjjkl — ^jikl — ^ijlk — ^jilk and ^*jkl — ^jikl ^ijlk ^jilk 5 

which is obtained by permuting the indices in pairs. 


Abbreviated Subscript Notation: 

Due to the symmetry of both the stress and strain tensors, there are at most 36 distinct elastic 
coefficients. Since many of the 81 coefficients are equivalent a reduced notation is possible. 
Using the abbreviated subscript notation developed by Voigt, we can rewrite the double indexed 
quantities as single indexed quantities as displayed in Table 1. 


Voigt Notation 

Ti = T„ 

£ 

II 

H 

III 

H 

T2 = T22 

t 5 = T 13 = T31 

T3 = T33 

T 6 = T 12 = T21 


s 4 = 2 S 23 = 2 S 32 

S2 = S22 

S5 = 2 S13 = 2 S31 

S3 = S33 

S6 = 2 S12 = 2 S21 


Table 1 


Following the common convention, we can write the elastic coefficients in reduced notation, taking 
into account the numerical factors, as 


C IJ ^ c ijkl 


and s D <=> s ijkl x 1 


1 for I and J = 1,2,3 

2 for I or J = 4,5,6 ► 
4 for I and J = 4,5,6 . 


(4) 


In the following discussion, quantities with lower case subscripts will refer to an 81 element cjj k i 
system while upper case subscripts, cn, will refer to the Voigt notation. Figure 5 illustrates the 
stress components acting on each face of a unit cube within a stressed body in the Voigt notation. 
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Figure 5: The stress components acting on each face of a unit cube within a stressed 
body (IJ Voigt system). 


The constitutive relation may be rewritten as 

6 

Ti = X c u S J (I = 1,2,3,4,5,6) 

J=i 

where cu represents the elastic stiffness coefficients. 

C 11 c 12 c 13 c 14 c 15 c 16 
C 21 c 22 c 23 c 24 c 25 c 26 
r I C 31 c 32 c 33 c 34 c 35 c 36 
C 41 c 42 c 43 c 44 c 45 c 46 
c 5 1 c 52 c 53 c 54 c 55 c 56 
_ c 61 c 62 c 63 c 64 c 65 c 66_ 


(5) 


If a localized strain-energy function exists then the elastic matrix that describes the solid is 
symmetric. Thus, the 36 elastic stiffness constants can be reduced to 21 distinct constants. A 
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material which has an elastic matrix of this form is called a triclinic system and represents a very 
general class of materials which exhibit no elastic symmetry. 

Reduction of Elastic Constants 

The number of distinct elastic coefficients needed to describe the elastic properties of a material 
can be reduced if the material exhibits elastic symmetries. For example, a large class of real-world 
materials exhibit orthotropic symmetry. These materials possess three mutually orthogonal planes 
of elastic symmetry and have at most 12 distinct elastic coefficients. Orthorhombic, tetragonal, 
hexagonal, cubic, and isotropic systems all exhibit orthotropic symmetry. We should note that the 
number of distinct coefficients will be reduced only if the coordinate axes chosen are aligned 
properly with the high-symmetry directions. Thus, an improper choice of coordinate axes will 
generally result in an increase in the apparent number of distinct elastic coefficients needed to 
describe the elastic properties of the material. In reality many of the coefficients will be linear 
combinations of the smallest set of coefficients obtained when the coordinate axes are chosen to be 
aligned with the principal axes. In Table 2 information is given for 7 linear elastic systems. 

The following discussion will demonstrate how a material possessing orthotropic symmetry 
reduces the number of elastic coefficients and produces an elastic matrix having a specific 
characteristic form when the material is aligned with respect to the principal axes of the system. 


Orthotropic System Derivation: 

The reduction of the number of distinct elastic coefficients can be illustrated by first 
considering a general rotation transformation a, 


lay] 


a ll a 12 a 13 
a 21 a 22 a 23 
a 31 a 32 a 33 J. 


(7) 


Application of this transformation to the stress and strain tensors yields the following relations 
between the rotated (primed) and nonrotated (unprimed) components of these tensors. 

T' = a T and S' = a S a^ (8) 


These relations can be written in matrix form as 


T' 


mn 


3 3 

^ y a mi Ty a n j 


i= 1 j = 1 


3 3 



i=l j = l 


(9) 


and 



Triclinic System: 21 coefficients (no 
elastic symmetries) 

C 11 c 12 c 13 c 14 c 15 c 16 
c 12 c 22 c 23 c 24 c 25 c 26 
j- -j c 13 c 23 c 33 c 34 c 35 c 36 
^ Cj 4 C 24 C 34 C 44 C 45 C 45 

c 15 c 25 c 35 c 45 c 55 c 56 
c 16 c 26 c 36 c 46 c 56 c 66 . 


Monoclinic System: 13 coefficients (only 
the xix 3 plane is a plane of elastic symmetry) 

Cn Cj2 C]3 0 c^ 0 

c 12 c 22 c 23 ® c 25 ® 

_ Cl2 C23 ° 33 ^ ° 35 ^ 

LcijJ = 0 0 0 C 44 0 C 46 

c 15 c 25 c 35 ® c 55 ® 

0 0 0 C 46 0 c 66 


Orthotropic Systems 

All systems below have 3 mutually orthogonal planes of elastic symmetry 


General Orthotropic Form: 12 

coefficients 

C 11 c 12 c 13 0 0 0 

C 21 c 22 c 23 0 0 0 

, _ C 31 ° 32 C33 ® ® ® 

lCljJ 0 0 0 C 44 0 0 

0 0 0 0 C<;<; 0 


0 0 0 0 C 55 0 

0 0 0 0 0 c 66 _ 


Tetragonal System: 6 coefficients 
(square packing) 

'c n c 12 c 12 0 0 0 

c 12 c 22 c 23 0 0 0 

_ Cl2 ° 23 C33 ® ® ® 

lCljJ = 0 0 0 C44 0 0 

0 0 0 0 C 55 0 

0 0 0 0 0 C 55 


Cubic System: 3 coefficients 

C 11 c 12 c 12 000 
c 12 C 11 c 12 0 0 0 
, _ c 12 c 12 C 11 0 ® 0 

LCtjJ 0 0 0 C44 0 0 

0 0 0 0 c 44 0 

0 0 0 0 0 C44 


Orthorhombic System: 9 coefficients 
(rectangular packing) 

C 11 c 12 c 13 0 0 0 

c 12 c 22 c 23 0 0 0 

C 13 c 23 C 33 0 0 0 

LCljJ “ 0 0 0 C 44 0 0 

0 0 0 0 C 55 0 

0 0 0 0 0 c 66 


Hexagonal System: 5 coefficients (six-fold 
axis aligned with xi axis) 

'cn c 12 c 12 0 0 0 

c 12 c 22 c 23 ® 0 0 

Ci 2 c 23 C 22 0 0 0 


[cn ] — 


c 22 “ c 23 
2 
0 


C55 0 
0 c 55 


Isotropic System: 2 coefficients 
Same form as the cubic system except 


c 12 — C 11 2 C 44 


Table 2: Linear elastic systems for which a localized strain-energy function exists. 
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S' 


op 


3 3 -> J 

- II a ok ^kl a pl II ^ok ^pl ^kl 


k= 1 1= 1 


k = 1 1=1 


( 10 ) 


Since we want the elastic stiffness coefficients to be invariant under the transformation, we make 
use of the constitutive relation and impose the condition that 

T’ = c’ : S’ => T = c : S’ (11) 


Equation (11) yields two sets of equations which when compared will determine which coefficients 
are to be eliminated. 


The following algorithm sketches the necessary steps to be performed in order to reduce the 
number of distinct coefficients needed to describe an orthotropic system. 


Step 1: 

Write out the primed (transformed) stress and strain elements in 3x3 form as a function of the 
unprimed variables. 


T' 


mn 


3 3 

^ , ‘hni Tjj 
i=l j=l 


3 3 

II ^mi a nj Tij 

i=l j=l 


( 12 ) 


S 


op 



k = l 1=1 


3 3 

II ‘kik ®pl ^kl 

k= 1 1=1 


(13) 


Step 2: 

Rewrite the sets of equations obtained in Step 1 using the Voigt notation. 


Step 3: 

Write the primed constitutive relation in component form. 

6 

T'i = X C H S J 

J= 1 

Replace the S'i's with equivalent Sj terms as given in Step 2. 

Replace the T'i's with equivalent Tj terms as given in Step 2. 

We now have the six primed equations totally in terms of the unprimed variables. 


( 14 ) 
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Step 4: 

Compare the set of equations from Step 3 with the original unprimed constitutive equations to 
determine which cu are to be set to zero. 



Planes of Elastic Symmetry: 

A plane of elastic symmetry exists at a point where the elastic constants have the same values 
for every pair of coordinate systems which are the reflected images of one another with respect to 
the plane. The axes of such coordinate systems are referred to as equivalent elastic directions. If 
the x lX3 plane is one of elastic symmetry, the constants cu are invariant under the coordinate 

transformation 

X 1 ! = Xj, x' 2 = - X 2 , x *3 = x 3 , (16) 

which maps a right-handed coordinate system into a left-handed system as illustrated in Figure 6. 



Figure 6: Geometry for illustration of an xjX 3 plane of elastic symmetry 

Let the xix 3 (or equivalently, x'ix’ 3 ) plane be a plane of elastic symmetry. Then 
Ti = cijSj andT’i = c n S’j ( J = 1,2,3, 4,5,6 ). The transformation matrix for mapping x 2 

to x' 2 is given by 
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[ajj] 


1 0 0 

0-10 


0 0 1 


(17) 


The stress and strain transformation laws yield the following relations. 


II 

t-H 

H 

T 11 , 

H 

<—> » 
to 

II 

-T 12, 

H 

II 

T 13 


T '21 = 

-t 21 . 

T '22 - 

T22> 

T ’23 “ 

-T 23 


T’ 3 l = 

T 31 . 

T '32 = 

-T 32 , 

T ’33 = 

T 33 

(18) 

S'n = 

S 11 , 

S'i2 = 

_ Si 2 , 

S'i3 = 

$13 


S*21 = 

_ S 2 i, 

$22 = 

$22’ 

S’23 = 

_ s 23 


S'31 = 

^31> 

S’32 = 

— ^32’ 

S'33 = 

$33 

(19) 


which imply 

T'i = T t 1 
S’i = sj 


for 1= 1,2, 3, 5 


and 



for 1 = 4,6 


( 20 ) 


Using the information in Equation (20) we are now able to write the primed variables totally in 
terms of the unprimed variables. 


Primed Constitutive Equations: 

T'j => Tj = c n Sj + c 12 S 2 + c 13 S 3 - c 14 S 4 + c 15 S 5 - c 16 S 6 

T ’ 2 => T 2 = c 21 Sj + c 22 S 2 + c 23 S 3 — c 24 S 4 + c 2 5 S 5 - c 2 g S 6 

T' 3 => T 3 = C 31 S 3 + c 32 S 2 + c 33 S 3 -c 34 S 4 + c 35 S 5 - c 36 S 6 

T ' 4 => T 4 = -c 41 Si - c 42 S 2 - c 43 S 3 + c 44 S 4 - c 45 S 5 + C 46 S 6 

T ' 5 => T 5 = C 51 S! + c 52 S 2 + c 53 S 3 - c 54 S 4 + C55 S 5 - c 56 S 6 

T’g =* T 6 = -c 61 Sj -c 62 S 2 -c 63 S 3 + c 64 S 4 -c 65 S 5 + c 66 S 6 (21) 

The relations given in Equation (21) must be equivalent to the unprimed stress/strain relations in 
the nonrotated coordinate system. 
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Unprimed Constitutive Equations: 

Ti = c n Sj + c 12 S 2 + c 13 S 3 + c 14 S 4 + c 15 S 5 + c 16 S 6 

^2 = C 21 Sj + C 2 2 s 2 + c 23 S 3 + c 24 S 4 + c 25 S 5 + c 26 S 6 

T3 = c 3i Sj + c 32 S 2 + c 33 S 3 + c 34 S 4 + c 35 S 5 + c 36 S 6 

T4 = C 41 Sj + c 42 S 2 + c 43 S 3 + C44 S 4 + C45 S 5 + c 46 S 6 

T5 = c 5i Si + c 52 S 2 + c 53 S 3 + C54 S 4 + c 55 S 5 + c 56 S 6 

^6 = c 6i Si + c 62 S 2 + c 63 S 3 + c 64 S 4 + c 65 S 5 + c 66 S 6 

Therefore, we obtain from Equations (21) and (22) that 

T 1 => C14 = C15 = 0 

T' 2 => c 24 = c 26 = 0 

T' 3 => c 34 = c 36 = 0 

T4 ^ C41 = c 42 = c 43 — C45 = 0 

T'5 => c 54 = c 56 = 0 

T'6 => C 61 = c 62 = c 63 ~ c 65 “ 0 


( 22 ) 


(23) 


and the elastic stiffness matrix reduces to 

C 11 c 12 C 13 0 C 15 0 
c 12 c 22 c 23 0 c 25 0 
r -j c 12 c 23 c 33 0 c 35 ® 

0 0 0 C44 0 c 46 

c 15 c 25 c 35 0 c 55 0 

0 0 0 c 46 0 c 66 _ ^ 24 ) 

Inspection of the elastic matrix representation for a material with one plane of elastic symmetry 

(xix 3 plane) yields the matrix form for a monoclinic system as displayed in Table 2. 


By repeating the procedure, this time requiring the xix 2 plane to be a plane of elastic 
symmetry, the following elastic stiffness coefficients can be shown to be zero. 
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T"! 


o 

Lti 

II 

o 

ON 

II 

O 

T" 2 

=> 

c 25 ~ c 26 - 0 

»TWI 

1 3 

=> 

c 35 = c 36 = 0 

, T’tt 

1 4 


c 45 = c 46 = 0 

T" 5 


C 51 = c 52 = c 53 = c 54 

T" 6 


C 61 = c 62 = c 63 = c 64 


(25) 


The elastic stiffness matrix for a medium having an x } x 3 and xix 2 elastic plane of symmetry 
reduces to 


C 11 c 12 C 13 0 0 0 

C 21 c 22 c 23 0 0 0 

C 31 c 32 c 33 0 0 0 

0 0 0 c 44 0 0 

0 0 0 0 c 55 0 

0 0 0 0 0 c 66 


(26) 


The elastic symmetry with respect to the (third) x 2 x 3 plane is identically satisfied by this 
matrix. That is, if we have two orthogonal planes of elastic symmetry with respect to principal 
axes, we must necessarily have the third. Inspection of Equation (26) shows that for an 
orthotropic material aligned along the principal axes, the elements of the upper-right and lower-left 
quadrants, and the off-diagonal elements of the lower-right quadrant are zero for the elastic matrix. 
The physical significance of these elements being zero can be understood by again returning to the 
constitutive relation. 

Inversion of the Constitutive Relation: 

To understand the physical role played by the elastic coefficients we must first transform the 
constitutive relation into a more convenient form. That is, we must obtain a form where we 
express the strain as a function of stress. This will allow us to systematically apply stresses along 
specific directions in order to ascertain which coefficients are involved in producing the strains and 
what types of strains result from a specific applied stress. We can invert the constitutive relation 
by multiplying each side of Equation (1) by the inverted stiffness matrix 

[c]" 1 : T = [c] _1 c : S = S => S = s : T (27) 

where [s] = [c] 1 is defined as the compliance coefficients. Equation (27) is written in matrix 
form as 
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Sjj = ^ ^ s ijkl T ki “ 1’2,3) 

k = 1 1=1 
6 

Si=I Si j t j (1=1,23,4,5,6) 


where for the most general case 


S 11 

Sl2 

s 13 

s 14 

S 1 5 

si6 

S 21 

s 22 

s 23 

s 24 

s 25 

s 26 

S 31 

s 32 

s 33 

s 34 

s 35 

s 36 

S 41 

s 42 

s 43 

S44 

s 45 

s 46 

S51 

s 52 

s 53 

s 54 

s 55 

s 56 

S 61 

s 62 

S 6 3 

s 64 

s 65 

s 66 . 


( 28 ) 


(29) 


Elucidating the Relationship Between Physical and Engineering Parameters 

We can obtain a better understanding of the physical significance each of the elements of the 
compliance matrix plays by investigating their dependence on the common engineering parameters 
for an orthotropic system. 


Young's Moduli: 

A Young’s modulus describes a physical phenomenon resulting from a single direction 
operation. That is, we apply a tensile/compressive normal stress along a given principal axis 
direction and observe the resulting strain along the same axis direction (see Figure 7). 

applied tensile/compressive normal stress 

Young's Modulus = resu j t j n g tensile/compressive normal strain (30) 

(both along the same principal axis direction) 


Consider, for example, the application of a tensile stress Ti (T n ) along the xi direction, with 


all other stresses zero. Then 

"sr 

5 2 

5 3 

5 4 

% 

_V 


S 11 s 12 S 1 3 0 0 

S 21 s 22 s 23 0 0 

S 31 s 32 s 33 0 0 

0 0 0 S 44 0 

0 0 0 0 S55 

00000 


0 

0 

0 

0 

0 

s 66 



'v 


0 


0 


0 


0 


_ 0 . 


(31) 
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Figure 7: Resulting strains produced by a tensile normal stress Tj (Tn). 


By inspection we have 

Si 

s 2 

S 3 


s n Ti 

(32. a) 

s 2 l T] 

(32. b) 

S 31 T 1 

(32. c) 


and S 4 = S 5 = S 6 = 0. The first of these equations allows us to express the Young's modulus 
in terms of the compliance matrix elements. 


Ei = 



By an analogous procedure we obtain 




(1= 1,2,3) 


(33) 


(34) 


As seen from Equation (34) the Young’s moduli are the reciprocals of the diagonal elements of the 
upper-left quadrant of the compliance matrix. 


Poisson Ratios: 

The Poisson ratios describe an operation which involves two orthogonal principal axes. We 
apply a stress along a given principal axis and observe the ratio of the resulting strain along a 
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principal axis perpendicular to the stress axis to the resulting strain along the stress axis (see 
Figure 7). 

resulting lateral strain “ 

resulting tensile/compressive normal strain . 

applied tensile/compressive normal stress) ( 35 ) 


Poisson Ratio = 

(due to an 


From Equation (32.b) we obtain 


V 21 



=> S 2 1 = 


S 21 T 1 
s il Ti 


v 2 i 

Ei 


f21 

S 11 


“ S 21 E i 


(36) 


A general relationship between the Poisson ratios and the compliance elements can be written as 



=> s jr = 


_ s ji T i 
s n T i 

Vji 

Er 


s n 

s n 


~ s ji e i 


( l,J = 1,2,3; i*J) 


(37) 


From Equation (37) we see that a specific Poisson ratio is obtained from a multiplication of the 
corresponding compliance element and the appropriate Young's modulus. 


Shear Moduli: 

A pure-shear strain results from the application of two diametrically opposing planar tangential 
stresses of equal magnitude acting on adjacent orthogonal surfaces to produce a pure distortion (see 
Figure 8). The direction of the applied forces are tangential to the surfaces to which they are 
applied. A shear strain is defined as 


Shear Modulus 


applied planar tangential stress 
resulting planar tangential strain , 


( 38 ) 


where the applied stress and the resulting strain are in the same plane. 
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Figure 8: Pure-shear strain in the X2X 3 plane 


Consider, for example, the application of planar tangential stresses T 4 (T23, T 32 ) in the x 2 x 3 
plane, with all other stresses zero. 

Sj 1 s 12 s 13 0 0 0 0 

5 2 S21 s 2 2 s 23 0 0 0 0 

5 3 s 3 j s 3 2 s 33 000 0 

s 4 = 0 0 0 S44 0 0 t 4 

55 0 0 0 0 S55 0 0 

5 6 _ _ 0 0 0 0 0 s 66 J L 0 J ( 39 ) 

Multiplication of the matrices yields the following stress/strain relation for a pure-shear strain in the 
X2X 3 plane, 

S 4 = s 44 T 4> ( 40 ) 

where Si = S2 = S 3 = S5 = S6 = 0 . From Equation ( 40 ) we see that the pure-shear modulus for 
the X2X 3 plane is equal to 

° 23 S 4 It 4 s 44 _ ( 41 ) 
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By an analogous procedure we see that the diagonal elements of the lower-right quadrant of the 
compliance matrix are related to the pure-shear strains in a reciprocal manner. 



T 4 , 

1 




1 


0 

to 

w 

111 

S7* T 4 

s 44 

:=> 

S44 


G 23 

(42) 


T 5 I 

1 




1 


G 13 = 

S5 It 5 

s 55 


s 55 


G 13 

(43) 


T 6 I 

1 




1 


G 12 = 

s 7 It * 

s 66 


s 66 


G 12 

(44) 

Combining the information obtained above we 

can rewrite the compliance matrix of an 

orthotropic system in terms of the engineering parameters as 
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Vl2 

Vi3 

0 

0 

0 



Ei 

E 2 

E 3 






V 21 

1 

v 23 
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0 

0 



Ei 

E 2 

E 3 






v 3 i 

v 32 

1 

0 

0 

0 


[s n ] - 

Ei 

E 2 

E 3 

1 





0 

0 

0 

G 23 

0 
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0 

0 

0 

0 

G 13 

- 0 

l 



0 

0 

0 

0 

0 

1 



_ 





^12 . 

(45) 


For orthotropic materials we have seen that many of the elastic coefficients are zero when the 
system is aligned with the principal axes. Some of these coefficients may not be zero if the 
material is not aligned. But since these coefficients will be linear combinations of the smallest set 
of distinct coefficients, the response of the material will also be a linear combination of the 
responses described above. 

For example, consider an orthotropic material prepared in the shape of a cube for which one 
side is aligned with the X2 principal axis while the two adjacent orthogonal surfaces are at an angle 
of 45° degrees with respect to the xj and X 3 principal axes of the system. Figure 9 illustrates the 
geometry of the material. By application of the Bond transformation matrix we obtain the 
following matrix for sjj. 



0 


s 'l 1 s 'l 2 s 'l 3 ® s 15 
s 2 1 s 22 s 23 0 s 25 ® 
s' 31 s ' 3 2 S' 33 0 S' 35 0 

ls ' IjJ " 0 0 0 s' 44 0 s' 46 

s 5 1 s< 52 s 53 ® s 55 ® 

0 0 0 s' 64 0 s' 66 J (46) 

Each of the s'u coefficients can be written as a linear combination of the smallest set of distinct 
compliance coefficients that describe the material when aligned with the principal axes. 





Figure 9: Application of a tensile normal stress T'i (T n) to an orthotropic material. 
The principal axes for the material is represented by the bold unprimed set of axes. 


Now consider the application of a tensile normal stress, T'i (T’n), along the x i axis with all 
the other stresses zero (see Figure 9). 


■s’r 


's'n 

S'l 2 

S'l3 

0 

S'l5 

0 

S2 


S'l2 

s ’22 

s> 23 

0 

S'25 

0 

S’ 3 


S’l 3 

s 23 

S' 3 3 

0 

S’3 5 

0 

S' 4 


0 

0 

0 

s '44 

0 

s 46 

S's 


s 'l5 

s 25 

s '35 

0 

s> 5 5 

0 

1 

CO 

CTs 

1 


_ 0 

0 

0 

s 46 

0 

s 66 



"TV 


0 


0 


0 


0 


. 0 . 


(48) 


Expanding Equation (48) the resultant strains can be explicitly written as 

S\ = s’ n T'j = |(s 11 + S 13 + S3 1 +S33 + 4S5 5 )T , 1 


S'n Tj 


( 49 ) 
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S'2 - s'21 T'j - ^ ( s 21 + s 23 ) T'j 


(50) 


S'3 - S 31 T'l - 7 ( S 11 + s 13 + S 31 + s 33 - 4 S55 ) T\ 


(51) 


S'5 - s 5i T'i - 7 ( Sji + s 13 - s 31 - s 33 ) T'j 


(52) 


and S '4 = S '6 = 0. Equations (49-52) can also be written in terms of the engineering moduli 
which more vividly illustrates that the response of the material is a combination of bulging, 
constricting, and shearing effects. 


S' - s' T - if -i-Iii 2 Ii + J_ + J_ It- 
1 " Sl,T ‘ “ 4U E, E, + E 3 + G 13 J T > 


(53) 


S- 2 = s' 21 T\ = - 7 


1 f V 21 v 23 

' 21 + -rr- I T'j 


2 l E 


1 


(54) 


S' = s’ T = ±- J_ It- 

3 31 1 HE, E 3 Ej + E 3 G 13 l Tl 


(55) 


S’s = s' 5 i T\ = 7 


1 v 13 v 31 1 


V E, 


T'i 


J 3 J 


(56) 


The final shape of the material is determined by the superposition of the contributions of all the 
nonzero strains (S'i, S' 2 , S' 3 , S' 5 ). Inspection of Equations (53-55) illustrate that the sole 
application of a tensile normal stress T'i along the x', axis produces normal tensile/compressive 
strain components along the x'j, x' 2 , and x ' 3 axes. The amount of the elongation or constriction is 
determined by the linear combination of the engineering moduli measured with respect to the 
principal axes. Equations (52) and (56) reveal that a shear strain also results from the application 
of a tensile normal stress under these conditions. 

As a second example, we will consider the case when all the stresses are zero except for planar 
tangential stresses T '5 (T'i 3 , T’ 3 i) in the x'ix ' 3 plane (see Figure 10). 
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Figure 10: Application of a planar tangential stress T ’5 (T’ 13 , T' 31 ) in an x'ix '3 
plane. The principal axes for the material is represented by the bold unprimed set of 


axes. 

The stress/strain relation is given by 
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(57) 

From Equation (47) we can write the primed compliance coefficients in terms of the unprimed 
coefficients. 


1 


S’l = S' 15 T’ 5 = T ( s ir s 13 + S 31 ~ s 33 ) T 5 


(58) 


S ' 2 = s '25 T ’ 5 = J ( S 21 - S 23 ) T ' 5 


(59) 
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S '3 = s '35 T 5 - j (s n - s 13 + s 31 - s 33 ) T ' 5 

S ’ 5 = s t 5 5 T '5 = ■j (sn - s 13 - s 31 + S 33 ) T 5 
Equations (58-61) can also be written in terms of the engineering moduli. 


s ’i - s 'i5 T'5 - T 


S ’2 = s' 25 T ' 5 = - - 


S’ 3 = s'35 T' 5 = T 


(60) 

(61) 
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(64) 
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Vl3 

1 

V31 
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1 > 

i TV 


Ei 

E 3 

Ei 

e 3 J 

! 1 5 

(65) 


S '5 - s '55 T' s - T I ~ + 


Again the superposition of the contributions from all the nonzero strains (S'i, S' 2 , S' 3 , S' 5 ) 
will determine the final shape of this material. Inspection of Equations (62-64) illustrate that the 
application of planar tangential stresses in the x'ix '3 plane (for this material) produces normal 
tensile/compressive strain components along the x'i, x' 2 , and x '3 axes. The amount of the 
elongation or constriction along a particular axis is determined by a linear combination of the 
engineering moduli measured with respect to the principal axes. Equation (65) illustrates explicitly 
how the shear strain contribution is dependent on the Poisson's ratios and Young's moduli for a 
sample, subjected to a shear stress, that is not aligned with the principal axes. 


Reciprocity Relation: 

Materials for which a symmetric strain-energy function exists exhibit a symmetric compliance 
matrix. 


Sjj — Sji (for I, J — 1,2,3, 4, 5, 6 ; I ^J) 

From Equation (36) for all stresses zero except Ti, we obtain 

S 2 | S 2 | _ v 21 

S21 = tT'Ti = s 7 e 7 T 1 “ Ej 


( 66 ) 


( 67 ) 


In a similar manner, if all the applied stresses are zero except T 2 , then 
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S 1 | S i I _ *12 

Sl2 “ T 2 It 2 S 2 E 2 It 2 E 2 

Therefore, if S21 = S12 then this implies 

V 21 _ v 12 

Ei _ E 2 


Thus, in general the equivalence of su = sji reveals that 

^ ^ (I,J = 1,2,3; I*J) 

b i Ej 


( 68 ) 


( 69 ) 


(70) 


Physical Interpretation of the Most General Compliance Matrix: 

The role of the compliance coefficients is to relate the applied stresses to the resulting strains 
for all possible combinations. Initially, we were required to have 81 compliance coefficients to 
relate uniquely the 9 stress components to the 9 strain components. In order to illustrate the role 
each compliance coefficient plays, we will systematically determine the type of stress to the type of 
strain each element uniquely relates. The most general form for the constitutive relation, in the 
Voigt notation, is 

S 1 s n s 12 s 13 s 14 s 15 s 16 Tj 

s 2 s 21 s 22 s 23 s 24 s 25 s 26 T 2 

s 3 s 31 s 32 s 33 s 34 s 35 s 36 T 3 

s 4 S 41 s 42 s 43 S44 S45 s 46 T 4 

S5 s 51 s 52 s 53 S54 s 55 s 56 T 5 

S 6 J L S 61 S 62 s 63 s 64 s 56 s 66 J L T 6 J . (71) 

Thus, the total resulting strain experienced by a material, due to the application of stresses, is the 
linear superposition of the contributions of Si, S2, S 3 , S4, S5, and S6- 

51 = s n T 1 + s 12 T 2 + s 13 T 3 + s 14 T 4 + s 15 T 5 + s 16 T 6 

52 = s 21 Tj+ s 22 T 2 + s 23 T 3 + s 24 T 4 + s-^ T 5 + s 26 T 6 

53 = S 31 Tj + s 32 T 2 + s 33 T 3 + s 34 T 4 + s 35 T 5+ s 36 T 6 

5 4 = s 41 Tj + s 42 T 2 + s 43 T 3 + S44T4+ S45T5+ s 46 T 6 

55 = S51 Tj + S5 2 T 2 + S5 3 T 3 + S54 T4 + S55 T5 + S56 T 6 

5 6 = s 61 T 1 + s 62 T 2 + s 63 T 3 + s 64 T 4 + s 65 T 5 + s 66 T 6 (72) 

Inspection of Equation (72) reveals that the compliance elements of the first column 
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s n for 1= 1,2,3.4,5,6 

relates the applied tensile/compressive normal stresses along the xi axis to the different types of 
strains. Table 3 summarizes the possible strain contributions resulting from the application of a 
tensile/compressive normal stress Tj (Tn) along the xj axis. 


Stress/Strain Contribution 

Resulting Physical Phenomenon 

S = Sjm Tn ^ ^1 — s n Ti 

Tensile/compressive strain along the xi axis 

S 12 = s 1211 Tn ^ ^6 = S 61 Ti 

Shear strain in the X 1 X 2 plane 

S 13 = s 1311 T lt => S 5 = s 51 Tj 

Shear strain in the X]X 3 plane 

S 2 I = S 21 11 Tji =* $6 = S 61 Tl 

Shear strain in the X 2 X 1 plane 

$22 = s 2211 ^11 ^ $2 = S 21 

Tensile/compressive strain along the X 2 axis 

$23 = s 2311 ^11 ^ $4 = S 41 Tl 

Shear strain in the X 2 X 3 plane 

B— B— 

Shear strain in the x 3 xi plane 

$32 — s 3211 Tn $4 _ S 41 Ti 

Shear strain in the x 3 X 2 plane 

$33 = s 33 1 1 Til ^ $3 = S 31 Ti 

Tensile/compressive strain along the x 3 axis 


Table 3: Contributions to the total resulting strain due to the compliance coefficients 
that involve Ti stresses. 


We see from Table 3 that the sn elements for I = 4,5,6 imply that the tensile/compressive normal 
stress T] would produce a shearing strain. For orthotropic systems these elements are zero when 
referenced to the principal axes for the material. Thus, a shearing strain cannot result from a 
tensile/compressive normal stress applied along a principal axis for orthotropic materials. 

From Equation (72) we see that the compliance elements of the sixth column 

s I6 for 1= 1,2,3,4,5,6 

relate an applied planar tangential stress in the xjX 2 plane to the different types of strains. Table 4 
summarizes the possible strain contributions resulting from the application of a planar tangential 
shearing stress T 6 (T 12 , T 21 ) in the X 1 X 2 plane. 















[r 


Stress/Strain Contributions 

Resulting Physical Phenomenon 

S u = s 1112 T 12 + s 1121 T 21 
=> Sj = s 16 T 6 

Tensile/compressive normal strain along the 
xj axis 

Si 2 = s 1212 T 12 + s 1221 T 21 
=> S 6 = s 66 T 6 

Shear strain in the x t x 2 plane 

s 13 = s 1312 T 12 + S 1321 T 2 i 
=> S 5 = s 56 T 6 

Shear strain in the X 1 X 3 plane 

S 2 j = s 2112 T 12 + s 2121 T 21 
=> s 6 = s 66 T 6 

Shear strain in the xix 2 plane 

s 22 = S 2212 T 12 + s 2221 T 21 
=> S 2 = s 26 T 6 

Tensile/compressive normal strain along the 
x 2 axis 

S 2 3 = S 2 3 12 T 12 + S 2 3 2 j T 21 

^ S 4 — s 4 g Tg 

Shear strain in the x 2 X 3 plane 

S 31 = s 3 1 12 ^12 + s 3121 ^21 
=> S 5 = S 56 T 6 

Shear strain in the X 1 X 3 plane 

S 32 = S 3212 T 12 + s 3221 T 21 
S 4 s 4 g T 6 

Shear strain in the x 2 X 3 plane 

S 33 = s 3312 T 12+ s 3321 T 21 

=> S 3 = s 36 T 6 

Tensile/compressive normal strain along the 
X 3 axis 


Table 4: Contributions to the total resulting strain due to the compliance coefficients 
that involve T 6 stresses. 


We see from Table 4 that S 66 relates a planar tangential shearing stress to a shear strain in the 
same plane. The s 46 and sse coefficients relate an applied shear stress in the xix 2 plane to shear 
strains in orthogonal planes, while the Si 6 , s 2 6 , and S 36 coefficients relate a shearing stress to a 
tensile/compressive normal strain along the principal axes. For orthotropic systems aligned with 
the principal axes S 46 = S 56 = s 2 6 = s 2 6 = S 36 = 0. Therefore, the strains involving these 
coefficients are not allowed. 

We can also view the compliance matrix in terms of quadrants (see Table 5), which groups 
together the coefficients by types of applied stresses and types of resulting strain contributions. 
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Tensile/Compressive 

Stress 

Planar Tangential 
Shear Stress 

Tensile/Compressive 
Strain Contribution 

S 11 s 12 s 13 
S 21 s 22 s 23 
S 31 s 32 s 33 

s 14 s 15 s 16 
s 24 s 25 s 26 
s 34 s 35 s 36 

Planar Tangential 
Shear Strain 
Contribution 

S 41 s 42 s 43 
S 51 s 52 s 53 
S 61 s 62 s 63 

s 44 s 45 s 46 

s 54 s 55 s 56 

s 64 s 65 s 66 


Table 5: The compliance coefficients grouped by types of applied stresses and the 
types of resulting strain contributions they uniquely relate. 


Figure 1 1 and Table 6 also illustrate the role each of the compliance coefficients plays in relating an 
applied stress to the resulting strain. 

As has been shown above, an orthotropic system has an elastic matrix which contains many 
coefficients equal to zero. If the material is not aligned with the principal axes, then these matrix 
elements may not be zero. These matrix elements are now linear combinations of the smallest set 
of distinct coefficients. (The Bond rotation formalism is a useful tool under these circumstances). 
Application of a tensile/compressive stress along the Cartesian axes of the misaligned system may 
produce both tensile/compressive strains and shearing strains. Thus, information on how the 
elements in the upper-right, lower-left and off-diagonal elements of the lower-right quadrants 
contribute to the resulting strains is important for materials whose macroscopic surfaces are not 
aligned with the principal axes for the material. 

In this Section, we have reviewed the formalism of linear elasticity in order to build a more 
solid foundation upon which an investigator can raise his/her physical intuition about the intrinsic 
nature of the physical properties of materials. We presented a physical interpretation for the role 
the elastic compliance coefficients play in uniquely relating applied stresses to resulting strains and 
the dependence these coefficients have on the traditional engineering moduli. Such information 
will be extremely useful in the design of advanced ultrasonic measurement systems. 
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Physical Interpretation of 
Elastic Coefficients 


Apply tensile/compressive normal stress; yield 
a tensile/compressive strain orthogonal to the 
applied stress axis 


- — • v.... ................. 


| Apply planar tangential shear stress; yield 
j tensile/compressive strain along an axis 
jorthogonal to the stress plane 


Apply tensile/compressive 
normal stress; yield a 
tensile/compressive strain 
along the axis of the stress 



s||l 

llll 

S 14 

S 15 

S 16 

-S> 


^23 

S 24 

S 25 

S 26 

S 31 

lllll 


S 34 

S 35 

S 36 

S 4 1 

S 42 

S 43 

^4^ 

^45 

S 46 

S 51 

S 52 

S 53 




S 61 

S 62 

S 63 

S 64 


|SI’ 


Apply planar 
tangential shear stress; 
yield planar shear 


Apply tensile/compressive normal stress; 
yield a shear strain in a plane orthogonal 
to the stress axis 


Apply planar tangential shear stress; yield 
planar shear strain in plane orthogonal to 
the stress plane 


S = s : T 


Figure 1 1 
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Applied Stress 
Involved 

Compliance 

Coefficient 

Valid Indices 

Resulting Strain 
Contribution 

Tensile/compressive 
normal stress along 
the xj axis 

SJJ 

J = 1,2,3 

Tensile/compressive 
strain along the xj 
axis, Axj/xj 

Tensile/compressive 
normal stress along 
the xj axis 

SlJ 

I, J = 1, 2, 3 
I * J 

Tensile/compressive 
strain along the xj 
axis, Axj/xj, an 
orthogonal axis 

Tensile/compressive 
normal stress along 
the xj axis 

SIJ 

I = 4, 5, 6 
J= 1,2,3 

Planar shear strain in 
the plane denoted by 
the I index, 

orthogonal to the xj 
axis 

Planar tangential 
shear stress in the J 
plane 

SJJ 

J = 4, 5, 6 

Planar shear strain in 
the same plane (J) as 
the applied stress 

Planar tangential 
shear stress in the J 
plane 

SIJ 

I, J = 4, 5, 6 
I* J 

Planar shear strain in 
the plane I, an 
orthogonal plane 

Planar tangential 
shear stress in the J 
plane 

SIJ 

1= 1,2,3 
J = 4, 5, 6 

Tensile/compressive 
strain along the X{ 
axis, orthogonal to 
the stress plane, 
Axj/xj 


Table 6: Physical interpretation of the elements of the elastic compliance matrix 
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